Overcritical Rotation of a Trapped Bose-Einstein Condensate 
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The rotational motion of an interacting Bose-Einstein condensate confined by a liarmonic trap is 
investigated by solving the hydrodynamic equations of superfluids, with the irrotationality constraint 
for the velocity field. We point out the occurrence of an overcritical branch where the system can 
rotate with angular velocity larger than the oscillator frequencies. We show that in the case of 
isotropic trapping the system exhibits a bifurcation from an axisymmetric to a triaxial configuration, 
as a consequence of the interatomic forces. The dynamical stability of the rotational motion with 
respect to the dipole and quadrupole oscillations is explicitly discussed. 
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An important peculiarity of harmonic trapping is the existence of a critical angular velocity, fixed by the oscillator 
frequencies, above which no system can rotate in conditions of thermal equilibrium. The main purpose of this work is 
to show that a trapped Bose-Einstein condensate at very low temperature can rotate at angular velocities higher than 
the oscillator frequencies in a regime of dynamical equilibrium. The occurrence of overcritical rotations is a rather 
well established feature in classical mechanics (see, for example, [|l],D) and is the result of the crucial role played by 
the Coriolis force. It is therefore interesting to understand the new features exhibited by rotating superfluids and in 
^ ■ particular the role played by Bose-Einstein condensation. The rotational behaviour of superfluids is, in fact, deeply 
, influenced by the constraint of irrotationality which makes it impossible, for such systems, to rotate in a rigid way. 
■ Spectacular consequences of irrotationality are the quenching of the moment of inertia with respect to the rigid value 
^ _ _ and the occurrence of quantized vortices |3|. Both these effects have been recently observed in dilute gases confined 
. in harmonic traps [||-||. 

' We start our analysis by considering a dilute Bose gas interacting with repulsive forces at zero temperature. For 
. large systems, where the Thomas-Fermi approximation applies, the equations of motion are well described by the 
so called hydrodynamic theory of superfluids [0J^. If the confining potential rotates with angular velocity Q it is 
, convenient to write these equations in the frame rotating with the angular velocity of the trap, where they take the 
form: 



rj^ . dp 

a : dt 



V(p(v- rj X r)) = (1) 



8: |.v.(^.^^.«#-v,.x,).„. 

- '—j In the above equations Vcxt(r) = M{uj'^x'^ + t^^y^ + t^z2^)/2 is the oscillator potential providing the external confine- 
ment. Notice that in the rotating frame T4xt(r) does not depend on time. Furthermore /i/oc(p) — gp is the chemical 
^ I potential of the uniform gas, where g = 47r fi^a/M is the coupling constant fixed by the s-wave scattering length a and 
- - ' V is the velocity field in the laboratory frame, expressed in terms of the coordinates in the rotating frame. It satisfies 
the irrotationality constraint. Eqs. ([^)-(^ can be applied also to a trapped Fermi superfluid where the expression 
for pioc{p) takes, of course, a different form. The stationary solutions in the rotating frame are obtained by imposing 
the conditions dp/dt = and dv/dt — 0. Let us look for solutions of the form |^ 

V = aV{xy) (3) 

for the velocity field, where a is a parameter that will be determined later. Choice (^) for the velocity field rules out 
the description of vortical configurations. Vortices cannot be in any case described by the hydrodynamic equations 
(|l|)-(^) since they require the use of more microscopic approaches, like Gross-Pitaevskii theory for the order parameter, 
accounting for the behaviour of the system at distances of the order of the healing length CoUq]. Since the critical 



frequency needed to generate a stable vortex becomes smaller and smaller as the number of atoms increases 1 1 1 1 , the 
solutions discussed in the present work correspond, in general, to metastable configurations. 

By substituting expression (^ into Eq. (^ one immediately finds that the resulting equilibrium density is given by 
the parabolic shape 



pir) = - 

g 



(4) 



1 



also in the presence of the rotation. Of course Eq. defines the density only in the region where p > 0. Elsewhere 
one should put p = 0. The new distribution is characterized by the effective oscillator frequencies 

Cbl ^ ul + - 2an (5) 
CjI ^ujl + a^ + 2an , (6) 

which fix the average square radii of the atomic cloud through the relationships 

CjUx') ^ Qliy^) ^ u;Uz') ^ ^ , (7) 



where the quantity 



^ 2 V hho ) 



2/5 

(8) 



is the chemical potential in the rotating frame and ensures the proper normalization of the density . In Eq. (H) 
we have defined Who = {oJx^y^jJzY^^ and aho — \J fi'/MCoho- The applicability of the Thomas-Fermi approximation, 
and hence of the hydrodynamic equations (^-(^, requires that the parameter 157Va/aho be much larger than unity. 
The rotation of the trap, providing a value of a different from zero, has the consequence of modifying the shape of 
the density profile, through the change of the effective frequencies Cjx and Cjy. For certain values of this effect can 
destabilize the system. Physically one should impose the conditions > 0, > to ensure the normalizability of 
the density. 

The equation of continuity (0), which at equilibrium takes the form (v — $7 x r) • Vp(r) ~ 0, yields the following 
expression for a 

/,r,2 _,r,2\ 

(9) 

in terms of and of the effective frequencies Lo^^ujy. From and (^) one finds that the expectation value of the 
angular momentum 

(L^) = M y"(r X v) ^ n(r) rfr = nQ (10) 

is always fixed by the irrotational value O — NM{{x'^ — y^))'^ / {x^ + y^) of the moment of inertia In terms of 
the effective frequencies lox and ujy the ratio between O and the classical rigid value 8rig = NM{x^ + 
simple expression 
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(11) 



Notice that both 8 and Qrig depend on the value of since the square radii {x^) and (y^) are modified by the rotation. 

Another useful quantity to calculate is the release energy E-c^i = -Ekin + -E-int giving the energy of the system after 
switching off the confining trap. This quantity can be extracted from time of flight measurements on the expanding 
cloud. In non rotating condensates it coincides with the interaction energy if one works in the Thomas-Fermi regime. 
In the presence of the rotation the kinetic energy Ekin = M J drn{r) cannot be instead neglected. By using the 
virial identity Q 2£'kin — 2i?ho + SEint — 0, where Eho is the expectation value of the oscillator potential, and noting 
that Eint — (2/7)Af/i, one finds the result 

z? . r, /, ,2 ,,,2\ 

(12) 

Let us now discuss the explicit behaviour of the stationary solutions of the hydrodynamic equations (0)-(^. By 
inserting expressions (||) and (|^) into Eq. (||) one finds the following third order equation for a jl^ : 

2a^ + a{ujl 40^) + n{ujl - w^) ^ q . (13) 
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Depending on the value of fl and of the deformation 



(14) 



of the trap, one can find either 1 or 3 solutions, derivable in analytic form. As already anticipated, the physical 
solutions should satisfy the additional requirements uj^ > and uj'y > 0, which ensure the normalizability of the 
density and rule out some of the solutions of (|l3|). The resulting phase diagram is reported in Fig. 1 where, in the 
plane — e, we show explicitly the regions characterized by 0,1,2 and 3 solutions (we assume hereafter e > 0). The 
solid curve, given by 

7 + —1-4^^ -] =0, 15 

divides the plane in two parts: on the left hand side Eq. ( p^ ) admits only one solution, on the right hand 
side it has three solutions. The dotted lines are the curves = tUy = Ux a/(1 -e)/(l + e) and fl = uj^. If 
e < 0.2 one can find (see Fig. 1) 3 stationary solutions, by properly choosing the value of fl |l4j. It is worth 
noticing that the phase diagram of Fig. 1 differs from the one derivable in the non-interacting Bose gas con- 
fined by the same harmonic trap. In this case the density profile in the rotating frame has the Gaussian shape 
p(r) = N {Mujho/T'fi')'^^^ Gxp [~M{CjxX^ + Cuyy^ + ujzz'^)/fi\. The renormalized frequencies still obey the Equations 
(^ and (^, but the relationship for a takes the different form a — —Vt{uOx — U0y)/{ujx +Cjy). In the non interacting 
case one finds that only one stationary solution is available for ^ < tOy and ^ > ujx, while no solution exists in the 
interval u)y < Q, < ojx- 

In Figs. 2 and 3 we show the stationary solutions of Eq. ( |l3| ) for a in two interesting cases: e = 0.1 where we 
predict the occurrence of a window with 3 stationary solutions, and e = 0.5 where one has a maximum of 2 stationary 
solutions satisfying the normalizability conditions a)^ > and > 0. The dashed-dotted line corresponds to the 
stationary solution of the non interacting gas. In both Figs. 2 and 3 one identifies two branches, hereafter called 
normal and overcritical branches. 

a. Normal branch: this branch starts at = 0. The linear dependence at small 17 is given by a = —fie. By 
increasing Q, the square radius {y^) increases and eventually diverges at O = ujy where Wy — > and the branch has 
its end |15|| . Also the angular momentum and the release energy diverge at 51 = LOy. Notice that when Wj, — > the 
moment of inertia takes the rigid value since {x^) <C {y^) (see Eq.(pT|)). 

b. Overcritical branch: this branch starts dXQ, — +oo where a behaves like a — (w^ — a>^)/4r2. It is worth noticing 
that in this limit both (D^ and ujy approach the value (w^ -I- u!y)/2 and therefore the shape of the density profile 
becomes symmetric despite the asymmetry of the confining trap. In the same limit the angular momentum tends to 
zero while the release energy approaches the finite value iJrei = (2/7)/ioo, where jloo is the chemical potential (||) with 
u;l — ujy = {luI + Ljy)/2. In the overcritical branch the deformation of the cloud takes a sign opposite to the one of 
the trap. This branch exhibits a back-bending at a value of f2 which is smaller than ujx, but can be higher or smaller 
than iUy, depending on whether the value of e is larger or smaller than 0.2 (see Figs. 2 and 3). In both cases this 
branch ends, after the back-bending, at the value = uix, where uix ^ and (a;^), (L^) and Ej-d diverge. 

It is also useful to discuss the instructive case e — s- corresponding to symmetric trapping in the x-y plane 
{ujx = LOy). In this case one finds a solution with a — for < lOx/V^ |16). For higher frequencies three solutions 
appear: the first one still corresponds to a non-rotating configuration (a = 0), while two solutions, given respectively 
by a = ±y2n^-wj, correspond to rotating deformed configurations. The existence of these solutions, which break 
the original symmetry of the Hamiltonian, is the analog of the bifurcation from the axisymmetric Maclaurin to the 
triaxial Jacobi ellipsoids for rotating classical fluids . It is worth pointing out that the existence of the bifurcation 
is the consequence of two-body interactions and is absent in the non-interacting Bose gas. When e is slightly different 
from zero (and positive) the two solutions are no longer degenerate, the one with a < having the lowest energy. 

The existence of stationary solutions in the rotating frame raises the important question of stability. Actually, 
one should distinguish between thermodynamic and dynamical instability. The former corresponds to the absence 
of thermodynamic equilibrium and its signature, at zero temperature, is given by the existence of excitations with 
negative energy p^ . The latter is instead associated with the decay of the initial configuration due to interaction 
effects and is in general revealed by the appearance of excitations with complex energy. Configurations characterized 
by thermodynamic instability are destabilized only in the presence of dissipative processes. 

Let us first discuss the stability with respect to the center of mass motion. In the presence of harmonic trapping 
the corresponding equations of motion, in the rotating frame, take the classical form (rotating Blackburn's pendulum, 
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|) and are not affected by interatomic forces. Their solutions obey the dispersion law 

1 



217^ ± 



y^(c.2_c.2)2+8f72(u;2+^2)J ^ (16) 



and are dynamically stable {uj"^ > 0) for fl < ojy and fl > ujx- So the requirement that the dipole oscillation be 
dynamically stable excludes the region between the dotted lines in Figs. 1, 2 and 3. Notice that this is the same 
region where the Schroedinger equation for the non-interacting Bose gas has no stationary solutions in the rotating 
frame. 

We have further explored the conditions of stability by studying the quadrupole oscillations of the condensate 
around the equilibrium configuration in the rotating frame. The calculation is derivable by linearizing the equations 
of motion (|l|) and (||), with the choice 

Sp{r) ^ ao + a.^x'^ + ayy^ + a^z^ + a^yXy (17) 
(5v(r) = V{axX^ + Uyy'^ + a-^z^ + a^yxy) , (18) 

for the fluctuations of the density and of the velocity field where the coefficients at and ai depend on time. The results 
of the analysis show that the window of dynamical instability for the quadrupole oscillations is different from the one 
of the dipole. In particular the normal branch is always dynamically stable, while the overcritical branch is stable only 
in its lower part since, after the back-bending, where da/dQ > 0, one of the quadrupole frequencies becomes purely 
imaginary {uj'^ < 0). We have investigated the stability of the quadrupole excitations also in the limiting case e = 0. 
In this case the upper and lower branches a = ±y^2fl'^ — lu'^ give rise to a vanishing value for one of the quadrupole 
frequencies, the others being always real. This vanishing solution corresponds to the rotation of the system in the 
X — y plane and reflects the rotational symmetry of the Hamiltonian. 

We stress again that the solutions discussed in this paper correspond, in general, to metastable configurations. For 
example it is well known that vortical configurations can become energetically favorable for relatively small values 
of Also excitations with higher multipolarities can become energetically favorable for some values of the angular 
velocity jl^. In this work we assume that at very low temperature, where collisions are rare, the solutions (^)-(|j) can 
survive, at least for useful time intervals, also in conditions of thermodynamic instability. 

Let us finally briefly discuss the experimental possibility of realizing the rotations described above. The normal 
branch, characterized by the huge increase of the size in the y direction and of the release energy when U, — s- a£j,, could 
be in principle generated by an adiabatic increase of the angular velocity, starting from a cold condensate [l|]. The 
transition from the normal to the overcritical branch cannot be instead realized in a continuous way and the only way 
to realize the exotic configurations of this branch is to engineer the proper conditions of dynamical equilibrium, by 
building up both the proper phase of the order parameter and the shape of the density profile. 

We are indebted to S. Vitale for introducing us to the problem of the overcritical rotations. Fruitful discussions 
with E. Cornell, J. Dalibard, L. Pitaevskii and G. Shlyapnikov are also acknowledged. This work has been supported 
by the Ministero della Ricerca Scientifica e Tecnologica (MURST). 
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FIG. 1. Phase diagram representing the stationary solutions of Eq. (^|) (see text). The dotted lines are the curves Q. = i 
and $1 = Ola; respectively. The full line is given by Eq. dl5b. 
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FIG. 2. Stationary solutions of Eq. ( p^ (solid lines) as a function of Q,, for e — 0.1. The dashed-dotted curves are the 
stationary solutions of the non-interacting Bose gas. The dotted straight lines correspond to Q, = uiy and Q, = lux respectively. 
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FIG. 3. As in Fig. 2, for e = 0.5. 
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